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Abstract 

The static and dynamic properties of the anisotropic XY-model (s = 1/2) on the inhomoge- 
neous periodic chain, composed of N cells with n different exchange interactions and magnetic 
moments, in a transverse field h, are determined exactly at arbitrary temperatures. The properties 
are obtained by introducing the Jordan- Wigner fermionization and by reducing the problem to a 
diagonalization of a finite matrix of nth order. The quantum transitions are determined exactly by 
analyzing, as a function of the field, the induced magnetization l/«X]m=i i^jm) ( 3 denotes 
the cell, m the site within the cell, a m the magnetic moment at site m within the cell) and the spon- 
taneous magnetization l/wX)m=l (^j,m>/ which is obtained from the correlations (^SJ m Sj +r m \ 
for large spin separations. These results, which are obtained for infinite chains, correspond to an 
extension of the ones obtained by Tong and Zhong(Physica B 304, 91 (2001)). The dynamic cor- 
relations, /Sj m (t)Sji m /(0)\, and the dynamic susceptibility, Xg 2 (^) ; are also obtained at arbitrary 
temperatures. Explicit results are presented in the limit T = 0, where the critical behaviour occurs, 
for the static susceptibility Xq Z (fy as a function of the transverse field h, and for the frequency 
dependency of dynamic susceptibility Xg Z (^)- 
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Keywords: anisotropic XY-model, quantum transition, inhomogeneous chain 



Electronic address: lindberg@fisica.ufc.br 



I. INTRODUCTION 



The one- dimensional XY-model introduced by Lieb, Schultz and Mattis [l| is still one of 
the few quantum many-body problems that can be solved exactly Although almost forty- 
five years have passed since the original solution has been proposed, this old model continues 
to provide new information on the quantum behaviour of mag netic systemsQ, 0, [J fl, fl. 
In particular, it has shed some light on the quantum phase transitions 7], since exact 
results can be obtained for most of its properties. The model has also been applied in the 
study of quantum entanglement, which plays an essential role in the quantum computation. 
Important results of this application can be found in the recent work by Amico et al. jlj and 
in the references therein. 

In this paper we will consider the anisotropic XY-model in a transverse field on the 
inhomogeneous periodic (closed) chain consisting of N cells composed of n sites, whose 
version on the open chain has been recently addressed by Feldman[9]. It has also been 
studied by Tong and Zhong lOj and Derzhko et al. ll], who have restricted their analysis 
to the study of the thermodynamic properties. To the best of our knowledege, the work by 
Tong and Liu 12j, which concerns the study of the zeros of the partition function, at T=0, 
and its relation to the quantum phase transitions contains the latest results on the model. 

The anisotropic model corresponds to an extension of the isotropic one recently studied 
by the anthorsQ. Its study on the inhotnogeneous periodic ehai„„a,™ponds to an 
extension of its version on the alternating superlattice|14j. 1151 . 



16 



17 



18 



19j . and, in this 



work, we will solve exactly the model by considering n different exchange constants and 
magnetic moments. The aim of this work is to present a comprehensive study of its static 
and dynamic quantum critical behaviour. 

In Section 2 we introduce the model and diagonalize its Hamiltonian. An explicit expres- 
sion is presented for n = 2, and the results compared to the known ones obtained by various 

19]. We also present the solution of the model for n = 8, which 
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authors 

has been obtained numerically. 

The induced magnetization M z and the isothermal susceptibility Xt '-> a ^ arbitrary tem- 
perature, are obtained in Section 3, and explicit expressions presented for T = 0. By ana- 
lyzing these quantities, the quantum critical behaviour is studied and the critical exponents 
associated to the quantum transitions are obtained. 
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FIG. 1: Unit cell of the inhomogeneous periodic chain 

The spontaneous magnetization M x is obtained in Section 4, from the two-spin correlation 
(Sf m Sf +rrr ^ , and the critical behaviour is also obtained, with high accuracy, numerically, 
by means of a non-linear regression of the data evaluated for finite values of r. This rather 
surprising result allows for the complete determination of M x as a function of the field, and 
to determine the multiple quantum transitions undergone by the model. 

The static and dynamic correlations (rf rf, r ) are presented in Section 5 and, the dynamic 
susceptibility x z *{uS), in Section 6. Finally, in Section 7, we summarize the main results of 
the paper. 



II. THE MODEL 



We consider the anisotropic XY- model (s = 1/2) on the inhomogeneous periodic chain 
with N cells, n sites per cell, and lattice parameter a, in a transverse field, whose unit 
cell is shown in Fig. 1, which corresponds to an extension of the isotropic model recentely 
considered [3]. The Hamiltonian is given by 



N 



1=1 {m=l 



n-1 



2_j (1 + 7m) Si m Si m+1 + J m (1 7 m ) Sf m Sf m+1 



m=l 



+ 



+ J n (1 + 7 n) S? !n S? +1>1 + J n (1 - 7 n) S^S^} 



(1) 



where the parameters J^ m are the exchange couplings between nearest-neighbour, \i m the 
magnetic moments, h the external field, and we have assumed periodic boundary conditions. 
If we introduce the ladder operators 



q± — q% _i_ qv 



in i 



(2) 
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and the generalized Jordan- Wigner transformation [2( 

l— 1 n m— 1 



K,m = eX P \™ C l,m' C l',rn> + ™ ^ 4™' C '> m ' f 



l'=l m'=l 
1 



m'=l 



(3) 



where ci.m and c| m are fermion annihilation and creation operators, we can write the Hamil- 
tonian as 21] 

H = H + P + + H~P-, (4) 



where 



N ( n 



H± = -E | (<=U- - 5) + E ft (■ 

Z=l l,m=l ^ ' m=l L 

1 1 JV_1 r j 

(4,m C !,m+l + C /,m+lQ,m) [ X Y ( C [n Q + 1 > 1 + 4+l,l C «,n) + 

' 1=1 L 

J m. f \ + 



m / f f , 

n I C l,m C l,m+l 1" C; m+1 Q iTn J 



JtrnTi', 



, ^»7n / t t 

+ (^ C i,n C i+l,l + Q+l,lQ,r 



± 



C jvn C l 1 + C 1,1 C ^,) 



and 



with P given by 



As it is well known 



I±P 



N 



P = exp [in c > 



(5) 



(6) 



(7) 



21 



22 



1=1 m=l 



231 ] . since the operator P commutes with the Hamiltonian, 



the eigenstates have definite parity, and P~(P + ) corresponds to a projector into a state of 
odd (even) parity. 

Introducing periodic and anti-periodic boundary conditions on c's for H~ and H + re- 



spectively, the wave- vectors in the Fourier transform [141 ] 



n 

\=^e^{-iqdl)a 
1 N 

= exp (igtZZ) ci irri 



(8) 



z=i 
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are given by q = for periodic condition and q + = ^j^p--, for anti-periodic condition, 



Nd 



with / = 0, ±1, ±N/2, and H and H + can be written in the form 



(9) 



where 



H 9 ± = -^2 Vmh (a\ ±m a q± , m -I)- 
m=l ^ ' 

n— 1 ^ 

- E ~2~ [ a 9 ± .™ a?± ' m+1 + a q ± ,m+i a 9 ± ,m + 7m ( fl g± , m a g± , m+ l + A^.m+l^.m) 



J, 



- JaJ ± n a g ± i iexp(-zdg ± ) + aj ± (1 a,±, n exp(zdg ± ) + 



+ 7m al± n al± i exp(-irfg ) + a^a^ exp(idq 



(10) 



Although if - and if + do not commute, it can be shown that in the thermodynamic limit 
all the static properties of the system can be obtained in terms of H~ or H + . However, 
even in this limit, some dynamic properties depend on H~ and H + [21, 22, 23 1. Since 
[H q ±, H_ q ±] 7^ 0, we consider the symmetrization 



H q ± = H q ±+H_ q ±, 



(11) 



and the Hamiltonian H can be written as 



(12) 



By making the identification q^ = q, it follows immediately that 



Hq, H q i 



0. Therefore, 



following the procedure introduced by Lieb, Schultz and Mattis we can diagonalize the 
Hamiltonian H ± by introducing the canonical transformation 



"q,m 



El Yq.k.m + Yq,k,m . Yq.k.m Yq,k,m + 
( ^ Vq,k + o VUk I 



a 



t 

q,m 



,m 0o,fc,m 



V-q,k 



(13) 



2 2 

where </> q) k,m and ip q ,k,m are the components of the eigenvectors, & qt k and Vf^fc, of the matrices 
(A, - B q ) (Aq + Bq) and (A, + B,) (A, - B,) , with A q and M q given by 
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2 



2 u 



exp 



2 



2 

Jz ■ . 

2 





V 4f exp 



/o 



Ji7i 
2 



-igd) 

Jl7l 
2 



J272 
2 





// h Jn ~ 1 



(14) 







J^rn exp 



J272 



J373 



J373 
2 



Jjh—\ T ri — 1 



exp (—iqd) 



Jn-17"-1 



(15) 



The corresponding eigenvalues are E^ k , the squares of the fermion energy levels, such that 



the Hamiltonian H , given in eq. ffTOl) . can be written in the diagonal form 



(16) 



where for H + (H ) the wave- vector q is identical to q + (q ), and the spectrum E q ^ of H q is 
determined from the determinantal equation 



det [(A q 



(A, 



(17) 



Since we are interested in the thermodynamic limit, we will consider from here on H~ 
only, in all calculations. 

Explicit expressions for the excitation spectrum can be determined for n equal to 2, 3 
and 4, and for n greater than 4, it is determined numerically from eq. fll7l) . Since we are 
considering the particle-hole representation, this implies that E q ^ 0. 

Explicitly, for n = 2, the excitation spectrum is given by 



El 



ai + a 2 + (Ai + A 2 ) cos 2q 



oil -a 2 + (Ai - A 2 ) cos2g] 2 + 4 [6f + S\ + 2<5i5 2 cos(2g)], 



(18) 
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where 



«i = (^) 2 + |(l+7i) 2 + |(l-7 2 ) 2 , 
« 2 = (M) 2 + #(l + 72) 2 + #(l-7i) 2 , 



A 1 = M(i + 7l )(i_ 72 ) ) 
A a =M(i_ 7l )(i + 72 ) | 

<*i = yC 1 - 7i)M + y (! + 7i)M, 

5 2 = y(l + 72)M + y(l-72)M- (19) 

For 71 = 72 = 0, the previous result reproduces the excitation spectrum obtained for the 
isotropic case 13] , provided we consider the particle representation which allows for negative 
energies. It also reproduces the result obtained by Siskens et al.[lf|. 

As already pointed out, although explicit expressions can be determined for n equal to 3 
and 4, they are not presented here since they are too cumbersome. 

Even in the case where we have identical magnetic moments, differently from the isotropic 
case, the effect of the the field does not correspond to a translation of the excitation spectrum 
for zero field, since the term of the field does not commute with the Hamiltonian. 

The critical fields which characterize the quantum transitions are determined from the 
excitation spectrum by imposing the condition E q = for q = and q = vr/dQ. The 
number of transitions is highly dependent on the anisotropy and varies from n to one. In 
particular, the limit of a single transition always occurs for the transverse Ising model which 
corresponds to the j's equal to one. 

The excitation spectrum for n = 8 and identical magnetic moments is presented in Figj2] 
for identical and different 7's, where we have considered the lattice spacing a = 1. As it can 
be seen, the energy of the modes is always positive and the effect of the spatial variation of 
the anisotropy is more pronounced in the low energy modes. Since there is no zero energy 
excitation, the model is not critical at the considered value of the transverse field. 
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FIG. 2: Excitation spectrum for n = 8, J\ = = J$ = 1.5, J2 = 3, J3 = 1.5, J4 = J5 = J7 = 2, 
h = 1, fi± = ... = fj,s = 1 and 71 = ... = 7s = 0.1 (continuous line), and for 71 = 0.1, 72 = 0.2, 
73 = 0.3, 74 = 0.4,75 = 0.5, 76 = 0.6,77 = 0.7, 78 = 0.8 (dashed line). 

III. THE INDUCED MAGNETIZATION AND ISOTHERMAL SUSCEPTIBILITY 

xf 



From eqs.Q, the operator Sf m can be written as 



S l, 



' 2 ^l,m-^l,m 1 



(20) 



with 



A^m — c\ + Cl )in , 



B, 



l,m — cj Cl tm , 



(21) 
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which can be expressed, from eqs.((8])and (fT3l) . as 



<i.l: 



By using eqs. (l20l and fl22l) . the local magnetization Mf m can be written as 

M tm = Mm <^ m > = -^E /im ^ fc - m ^^» ~ 2ra 9.fc) ' (23) 

q,k 

where the fermion occupation number n q ^ is given by 

n 9,k = - 7 , (24) 

1 + 

where E q u = f3E q ^ = E^k/ksT, and the calculation has been done by considering H = H 



21 



22 



231 ] . since we are interested in the thermodynamic limit. 



As in the isotropic model 



131 ]. we define an average cell magnetization operator in the z 



direction, if, as 

1 - 

r ^ ^ ^ a2 fe^mi (25) 

m=l 

then the induced magnetization per site, M z , is equal to (rf) and, from eq. (l23l . can be 
written in the form 

M z = (rf ) = ]T // m ^ )feim 0; fe , m tanh ( ) , (26) 

q,k,m 



which at T = becomes 



M z 



(27) 



2niV 

q,k,m 



The isothermal susceptibility can be obtained from eqs. (1261) and (|27|) by means of the ex- 
pression 

1 dM z 

Xt = (28) 
n oh 

which has to be evaluated numerically. 

At T = 0, the isothermal susceptibility Xt diverges at the quantum critical points induced 
by the field. This result can be seen in Figj3] where M z and Xt are presented as functions 
of the field h, at T = 0, for a chain with n = 8, identical /x's, 7 = 0.1 and 0.15. The 
magnetization, differently from the isotropic model, does not present plateaus but it does 
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FIG. 3: (a) Magnetization M z and (b) isothermal susceptibility T = 0, as functions of the 

field for n = 8, fix = ... = us = 1, J\ = Jq = Js = 1.5, J2 = 3, J3 = 1.5, J4 = J5 = J7 = 2, 
7 = 0.1(dashed line) and for 7 = 0.1 5 (continuous line). 



present inflexion points which induce the divergences of the susceptibility at the quantum 
phase transitions. As expected, the results also show that for n even there is the tendency 
of formation of a zero magnetization plateau in the limit 7 — > 0, which is not present for n 
odd as shown in FigHJfor n — 5. 

At T = 0, on the critical region, we have verified numerically that the isothermal suscep- 
tibility presents the behavior 

X¥~ln\h-h c \ (0< 7m <l), (29) 

as it can be seen in the results presented in Figj5l for a chain with n = 4, and identical //s 
and j's. 

The transitions belong to the same universality class as the one in the homogeneous 
model 24|. For T > all these transitions are suppressed by the thermal fluctuations. 
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FIG. 4: (a) Magnetization M z and (b) isothermal susceptibility T = 0, as functions of the 

field for n = 5, [i\ = ... = ^5 = 1, J\ = 1, J2 = 1.5, J3 = 2, J4 = 2.5 J5 = 3, 7 = 0.1(dashed line) 
and for 7 = 0.3 (continuous line). 



As already pointed out, the critical fields are associated with the zero-energy modes with 
q = and q = ir/d, and the number of critical points can vary from one in the transverse 
Ising model limit, for equal j's = 1, to n, provided the //s are positive. On the other hand, 
if we allow the fi's to assume negative values, we can have more than one phase transition 
even in the transverse Ising model limit, as has been shown by Derzhko et al. 11]. 

Then, the critical fields can be obtained numerically from the excitation spectrum and 
the results are shown in Fig-d^k) for a chain with n = 4, and identical /x's and j's. As can be 
verified, the number of transitions depends on 7 and its signature is present in the behaviour 
of the spontaneous and induced magnetizations. From these results, we can also conclude 
that for n even we have two critical j's, which can induce quantum transitions at zero field, 
whereas for n odd we have just a single critical 7, equal to zero, as in the uniform model 
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FIG. 5: (a) Magnetization M z and (b) isothermal susceptibility Xt ■> at T = and n = 4, \i\ = 
^2=^3 = ^4 = 1, 7i = 72 = 73 = 74 = 0.1, as functions of the uniform field for J\ = 1, Ji = 1.5, 
J3 = 1.75, J4 = 2.0. The logarithmic singular behaviour of Xt ^ ^ s shown in the inset in (b) where 
the horizontal scale is decimal logarithmic. 



251 ] . This result is shown in Fig.([7j) and, differently from the isotropic model, it constitutes 



the main effect of the cell size on the quantum critical behaviour of the model. 



IV. THE MAGNETIZATION M : 



Since the Hamiltonian is invariant under the transformation Sf m — > —Sf mJ which means 
that the < Sf m > is equal to zero, we cannot calculate the spontaneous magnetization 
directly by using the Hamiltonian shown in eq.([T]). On the other hand, the Hamiltonian 
after the introduction of an external field along x direction, which would eliminate the 
mentioned symmetry property, is no more exactly soluble. Therefore, by using the original 
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Hamiltonian, the spontaneous local magnetization, M^j, can be determined, from the static 
correlation function (Sf m Sf +r m ) , by means of the well known decomposition 26] 



lim (Si Sf + ) — \Sf m )(Si + ) . 
Then we can write immediately the local spontaneous magnetization as 



(30) 



M x m = lim J(Sf tm S? +r J, (31) 
and from this we obtain the average spontaneous magnetization per cell, M x , given by 

1 n 

MX = -J2 M ™- ( 32 ) 

m=l 

Therefore, by using the eqs.([2]), ([3]) and (l2~Tj) . the correlation (Sf m Sf +rm } can be written 



as 



< ^f,m^f+r,m ^ ^ ^,m^,m+l^,m+l^,m+2-^,m+2 • • • -Az+r.m— l-^+r,m-l^+r,m > • (33) 



Following [l| we can write (Sf m Sf +rrn ) , by using using the Wick's theorem, in terms of the 
determinant 



^ (Bi im Ai ;Tn+ i) 

(Bl t m+lAi tm+ i) 
(-B/,m+2^Z,m+l) 



(-Bz,m+2^,m+2) 



,m+l 



(B^ m A^ n ) 
(Bi 

(-Bz,m+2^4/,n) 

(Bi in Ai in ) 



(Bi >m Ay, m ,) = ± J2 e ~ t9d{l '' l) ^^rn^n^ - 1). 



(Bi, m Ai +r ^ m ) ^ 
(Si 

\-E?i,m+2^4z+r,m) 
\Bl n Ai+ r>m ) 



\ {Bl+ r) m-lAi yTn+ i) (Bi +r ^ m ^iAi tm+ 2) ■ ■ ■ {Bi +r , m ^iAi +r)n ) . . . (Bi +r ^ n Ai +ryTl ) J 

(34) 

where the static contractions (B^ m Aii . m /) , at arbitrary T, are given by 



(35) 



Therefore, in order to obtain the average spontaneous magnetization per cell, M x , we have 
to determine the asymptotic behaviour of the determinant shown in eq. (134"]) . For finite 
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temperature, it can be shown that the asymptotic behaviour of the determinant is zero, which 
corresponds to the limit r — > oo, and, as expected, there is no spontaneous magnetization 
at T + 0. 

On the other hand, at T = 0, the asymptotic behaviour of the determinant is different 
from zero for values of the field where the ground state is ordered. This asymptotic behaviour 
of the determinant can be estimated numerically by using the e— algorithm 2^], and, by 
using this method, we have obtained M x for a lattice with n = 4, identical y!s and j's. In 
order to obtain the estimate of the asymptotic value of the determinant, we have considered 
a numerical series constructed by varying r from 56 to 80, and the results are shown in 
Fig-dSb). The behaviour of the magnetization close to the transition points has been obtained 
by adjusting the reliable numerical results in the critical region to the scaling function 

M x {h) = M \h-h c f. (36) 

By considering a non-linear regression we have been able to obtain the adjustable parameters 
of the scaling function, namely, the critical field h c , the critical exponent (3 and the amplitude 
Mq. These results are shown in Table I. 

As it can be seen in Table I, the comparison of the numerical results, obtained from the 
non-linear regression, with the exact known ones show that they are extremely precise. This 
rather surprising result means that the quantum critical region can be precisely described 
numerically by eq. fl36|) . 



V. STATIC AND DYNAMIC CORRELATIONS {t^t^) 



The correlation function < Sf (t)Sf +r (0) >, in the thermodynamic limit, can be obtained 



from the expression 



21 



22 



23] 



/qz mqz , n xv _ Tr [ex P (-(3H-)ex P (tH-t)S^ex P (^H-t)S! +r ^] 
{b l;m it)b l+r)m ,{0)) - ______ . {67) 

Therefore, the dynamic correlation between the effective spins in the field direction, 



1 n 

(ri(t)r l %(0)) = - Mm'(W^(0)) 



(38) 



m,m'=l 
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FIG. 6: (a) Critical lines, at T = 0, for n = 4, /ii = fi2 = ^3 = \i± = 1 and for Ji = 1, J2 = 1.5, 
J3 = 1.75, J4 = 2. (b) Spontaneous magnetization M x as a function of the uniform field, for 
different 7, n = 4, /xi = Hi = 113 = (14 = 1 and Ji = 1, J2 = 1.5, J3 = 1.75, J4 = 2. The continuous 
line represents the extrapolation procedure. 

by using eqs. (j3]) and ( !2lT) . and Wick's theorem, can be expressed as 

j n n 

(?f (*K+r(°)) = ^ tWn* {( A l,m(t)Bi >m (t)) (A l+rtm/ (0)B l+r>m ,(Q)) - 

m=l m'=l 

— (A,m(*)A+r,m'(0)) (-B J)m (t).Bj +r!m /(0)) + 

+ (^ iffl (t)B (+r , m( (0)) (^^(^^/(O))} , (39) 
where the dynamic contractions, obtained from eq. ( 122ft . are given by 
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FIG. 7: Critical lines, at T = 0, for n = 5, /xi = fi2 = fJ-3 = A*4 = = 1 and for Ji = 1, J2 = 1.5, 
J3 = 2, J4 = 2.5, J5 = 3, where the vertical scale is decimal logarithmic. 





(0)) = 


1 

N 


q,k 




[e iE i' kt n q , k + e~ iE «> kt (1 


- w?,fc)] , 


\Bl tm (t)B l+r ^ m i 


(0)) = 


1 

N 


q,k 


-iqdri i* 

H J q,k,mH J q,k,m' 


[- e iE «- kt n q , k - e~ iE "- kt 


(! -w?,fc)] , 


(Al t m(t)Bi + r tm i 


(0)) = 


1 

N 


q,k 


t9dr ^q,k,m(i)* q: k,m' 


[- e iE *' kt n q , k + e~ iE "' kt 


(!-«?,*)] , 


(B^ m (t)Ai +r ^ m i 


(0)) = 


1 

N 


q,k 


tqd " '" ( l ) q,k,m^'q > k,rn' 


[e iEqM n q , k ~ e~ iE "- kt (1 
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Then we can write the dynamic correlation (t? (t)rf +r (0)) in the form 

(r l z (t)r l z + M) = (ri z ) 2 + 

q,k,m g',k',m' 

x [e^+^'^n^n^k, + e -^ q , k +E q , ik ,)t (i _ nq k) (i _ n ^ )fc/ ) + 
+ e*^'-*-^.*'^,,* (1 - n^) + e - i ( J5 «.*- E v,*') t (1 - n q , k ) n^] - 

q,k,m q' ,k',m' 

x [e i{E ^ +E "'- k,)t n q ,kn ql , kl + e -^ q , k +E q/ <kl )t (1 _ n ^ (1 _ ^ # ^ _ 
_ e i(s„ fc -^, t ,)t Nifc (x _ Uq , k ,) _ e -i(B Blfc -^, t ,)* (i _ Uq k) Uq , ^ . 

For t = we obtain the static correlation, which is given by 



<^(0K +r (0)) = (rf) 2 + 

q,k,m q',k' ,m' 
q,k,m q' ,k',m' 

x (1 - 2n 9;jfc ) (1 - 2n,/ )fc /) . 
At T = 0, the dynamic and static correlation are given respectively by 



r l z (t)r l U0)) = (r i r + 



An 2 N 2 

q,k,m q',k',m' 



- 4n 2 jv2 Y Y ^^' e_ * ( ^ 9 ' )dr ^,fc^^,fc,m'^'/ , ,fc' I ™'^',fc',m X 
q,k,m q',k',m' 

x [- e -i(B,, fc +^ ifc ,)t] 
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FIG. 8: Disorder lines (continuous line) and critical (dashed line), at T = 0, for n = 4, = fj>2 = 
Hz = M4 = 1 an d for J i = 1, J 2 = 1-5, J3 = 1.75, J4 = 2, where the vertical scale is decimal 
logarithmic. The asymptotic behaviour alternates between monotonic and oscillatory as we go 
from region I to region V. 



and 



<^(0K +r (0)) = (rf) 2 + 

q,k,m q' ,k',m' 
q,k,m q' ,k',m' 



As in the homogenous model 



25], the asymptotic behaviour of the static correlation 
(rf (0)r^ r (0)), at T = 0, as a function of r, can be oscillatory or monotonic depending on 
the values of the interaction parameters. For arbitrary n, these regions are separated by 
hypersurfaces in the parameter space and are called disordered surfaces, which collapse into 
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FIG. 9: Static correlation function in the field direction, at T = 0, for n = 4, fix = n% = fj,s = fj,^ = 1 
and for J\ = 1, J2 = 1.5, J3 = 1.75, J4 = 2, as function of r (distance between cells) for values 
of the field belonging to different regions shown in Fig. 8, where the vertical scale is decimal 
logarithmic. 



disorder lines for identical 7's. For n = 4, these disorder lines and the critical lines are shown 
in Fig. ([SD - Independently of the number of sites in unit cell, in the plane 7 X h, the disorder 
lines are given by equation 

^ + 7^ = 4, (45) 
d. Along these curves the anisotropic model is equivalent to an 



where hj C is a critical fie 
effective isotropic model [25j. For n = 2, this equivalence is shown in Appendix A, and for 
n > 2 it has been verified numerically. It is also shown, that, along this curve, the effective 
isotropic model is in a disordered state, which corresponds to the magnetization plateaus 
for the case where we have identical fi's. 

The asymptotic behaviour of the static correlation function, also for n = 4 and at T = 0, 
as a function of the distance r between cells, is presented in Fig.® for various values of h 
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FIG. 10: Asymptotic behaviour of the direct static correlation function in the field direction, at 
T = 0, iov n = 4, fj,i = fj,2 = fJ,3 = m = 1 and for J\ = 1, = 1.5, J3 = 1.75, J4 = 2, as function of 
r (distance between cells) for values of the field belonging to the critical lines (a), where the scales 
are decimal logarithmic, and disorder lines (b), where the vertical scale is also decimal logarithmic. 



and 7 in the different regions shown in Fig.®. As can be seen, the behaviour, as expected 
alternates between oscillatory and monotonic. 

At the critical point, as in the homogeneous case, the direct static correlation behaves 



asymptotically as r 



251 ] . This behaviour is shown in Fig.(fT0k). for n = 4, and different 



critical fields. However, at the disorder line, differently from the homogeneous model where 
the direct static correlation is zero, it behaves asymptotically as A _r , where A is a function 
of the field and of the interaction parameters, as shown in Fig.fllOb). for identical 7's and 
different values of the field. 

The real and imaginary parts of the dynamic correlation for a chain with n = 4 and 
uniform fi's, at T = 0, are presented in Fig. ( TTTT) for values of the field greater and equal to 
the critical field. As in the isotropic case, apart from the asymptotic behaviour, no noticeable 
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difference is observed in the dynamic correlation at different values of the field. 



VI. DYNAMIC SUSCEPTIBILITY 



From the dynamic correlation (rf (t)T* + r (0)) , eqs. (|4ip . we can obtain the time Fourier 
transform of dynamic correlation { T i T i+ r ) ^ from the equation 



oc 



W^> w = ^ / W(*)^Ar(0)> e**df, (46) 



-oo 



which is given by 



+ VmV m >e^ {q ~* )dr i>q,k,m^lk,m'&q',k',m4q>, k >,m X 



An 2 N 2 

q,k,m q 1 ,k',m' 



X [5 (UJ + E Qjk + £y,fc') n q ^nq',k' + H 1 ^ ~ E q,k ~ E q > >k >) (1 - 7lg )fc ) (1 - Tty,*/) + 

+ 5 (u; + E 9jfc - ^.fc/) n, )Jfc (1 - ry jfc >) + 5 (u - E q<k + i^,*./) (1 - n q<k ) n q , )k ,} - 

~ An 2 N 2 2.^1 ^rn^m'e YqArn4 > q,k,m'' l P<l',k',m' ( Pq',k',m X 

q,k,m q',k',m' 

x [5 (u + E q>k + Eq, ik >) n q> kn q ' t k' +5(u- E q>k - E ql>k ,) (1 - n q , k ) (1 - n q > >k ') - 
-6(u + Eq tk - E qlj k>) riq tk (1 - n q > ;k >) -6(u- E q>k + E q > >k >) (1 - n ff)fc ) ry jjfc /] . (47) 

By introducing the spatial Fourier transform 

W^^EW^)/**'. (48) 
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in the previous expression, we can write immediately 



W^> w = N5 q>0 6(u) (rn 2 + 

+ a 2 AT f I m^m'^q',k,m4'q^k,m' ( f ) q'-q,k',m'4 ) q'-q,k',rn X 



q' ,k,m,k' ,m' 



x [5(u + E q ' ik + E q i_ q ^ k i) n q ' jk n q '- q> k' + 5 (uj — E q ^ k — E q /_ q ^ k >) (1 — n q ^ k ) (1 — 



? -q, 



+ S (uj + Ey^ — E q '- q> k') n q > t k (1 — n q '- q> k') + 5 (u — E q > >k + E q i- q>k >) (1 — n g '_, 

_ 4fl 2 N / im ^ m '^9^ fc ." ^< ^g',fc,m'' ^ / ; g'-g,fc',m'0g'- (? ,fc',m X 

x [5 (u; + + E q i^ q ^ k i) n q f t kTi q '- q> k' + 5 (uj — E q i jk — E q - q / tk >) (1 — ?V,k) (1 — n. g -y, 
— 5 (a? + £y,fe — Eq/- q> k') n q > tk (1 — n q i- qtk >) — 5 (uj — .Ey^ + E q i^ qtk i) (1 — n^,*.) n 9 '_, 

(49) 



The dynamic susceptibility y^Huj) can be obtained by using the expression 28] 



UJ — UJ' 



(50) 



and from this we obtain 



4n 2 iV 



2 t» j ^ j f^rn^m' (Pq' ,k,rn,q' —q,k' ,m' Q q' ,k,in,q' —q,k' ,m') ^ 



X 



q' ,k,m,k' ,m' 
1 



W + E q ' jk + E q '- Qjk > UJ — E q > )k — E q >_q )k / 
~A 2 AT ^ . l^ml^m! (Pq' ,k,m,q' —q,k' ,m' Qq' ,k,m,q' —q,k' ,m' 



x [n q > jk + n q >- qik ' - 1] 



x 



x 



q' ,k,m,k' ,m' 
1 



u + E q i jk — E q i__ q k i uj — E a i k + E a i_ a k i 



q',k ~r ^q'-q,k' 



where 



Pq',k,m,q'-q,k',m' ~ VV ,fc,"i^V,fc,m' < A?' '-g,fc' \m'4 ) q> '-q„k> ' ,mi 
Qq',k,m,q'-q,k',m' = q 1 .k.m^ q\k,ra'' X l ) q 1 -1,k' ,ra' 4> q> -q^k' .rai 

and we have used the identity 

1 - n 9)fe = e pE "' k n qjk . 



(51) 

(52) 
(53) 

(54) 



It should be noted that the previous result reduces to the known one for the isotropic 



model 
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The static susceptibility Xq Z (fy is obtained by making u = in eq. (l5T]) . and is given by 

X,q (0) ~Z 2 t> j ^ . /^m/^m' [Pq' ,k,in,q' —q,k' ,m' Qq' ,k,m,q'—q,k' ,m') X 



2n 2 iV 

q' ,k,m,k' ,m' 

x ( — ^ ) x + n ql -q yk ' - 1] 



Eq',k + Pq'-q.k 1 
7~j\T2 ^ ^ i^Pq' ,k,m,q'—q,k' ,m' Qq' ,k,m,q'~q,k' ,m') X 



2n 2 iV 2 

q' ,k,m,k' ,m' 



X 



X K'.fe ~~ TV-?,*'] • ( 55 ) 



.-Ey.fc — Pq>-q,k' , 

The isothermal susceptibility can also be obtained by using the dynamic correlation given 



in eq. (j4T]) . and can be written as [29] 

*T=J2j - (rf) 2 ) dX, (56) 

r o 

where, as defined previously, (3 = l/ksT. From this it can be shown that the isothermal 
susceptibility in the field direction is equal to the uniform static one XiT(O)- 

At T = 0, the dynamic and static susceptibilities given in eqs. (l5"Tj) and (1551) can be 
explicitly written as 



(^°0 A.rfl Af ^ ^ t^mf^m' i,Pq' ,k,rn,q'—q,k' ,m' Q q' ,k,m,q' ~q,k' ,m') X 

(58) 



An 2 N 

q' ,k,m,k' ,m' 



1 



and 



LO + E q > t k + E>q'-q,k' W — -Eg', A: ~ E>q'-q,k' 



X,q (0) q 2 a 7" ^ ^ ^mf^m' i^Pq' ,k,m,q'—q,k' ,m' Q q' ,k,m,q'—q,k' ,m') X 



2n 2 N , 

g ,k,m,k ,m' 

E + £ (59) 
The static susceptibility Xg Z (0), at T = and for n = 4, as a function of the field is 
presented in Fig. ffT2]) . For q = 0, as expected, it diverges at the critical fields and it tends 
to zero as h — > oo, which corresponds to the saturation of the model. However, for q ^ 
differently from the isotropic model, no divergence is present, even for a wave- vector at the 
zone boundary. 

The real and imaginary parts of x^(u;) are obtained by considering x^i^* ~ * e ) i n t ne 
limit e — * in eqs. floT]) and fl55l) . The results, at T = and for n = 4, are shown in 
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FIG. 11: (a) The real and (b) imaginary parts of the dynamic correlation function (r^(i)r^ 1 (0)) , 
at T = for n = 4, /xi = /X2 = /U3 = /X4 = 1 and J\ = 1, J2 = 1.5, J3 = 1.75, J4 = 2, as function of 
time t. 



Fig. (fT51) for different wave- vectors, as functions of u. As expected, since the Hamiltonian 
of the model preserves the symmetry of the spin interactions of the homogeneous one, the 
divergences in the real part, for any wave- vector, correspond to square- root singularities in 
the imaginary part. 



VII. CONCLUSIONS 



In this work we have considered the anisotropic XY-model on the inhomogeneous periodic 
chain with N cells and n sites per cell. The model has been exactly solved, at arbitrary tem- 
perature, for the general case where we have n different exchange constants and n different 
magnetic moments. The number of branches of the excitation spectrum is equal to n, and 
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FIG. 12: Static susceptibility in the field direction, Xg 2 (0)i at T = 0, as a function of the field for 
n = 4, Hi = /j,2 = H3 = m = 1, J\ = 1, J2 = 1-5, J3 = 1.75, J4 = 2 , 7 = 0.1 and different values of 



analytical results can be found for n ^ 4. For n = 2, explicit expressions are presented for 
the excitation spectrum, and, for n — 8, it is obtained numerically. 

The induced magnetization M z and the isothermal susceptibility Xt are a l so given by 
explicit expressions at arbitrary temperature. At T = 0, where the quantum transitions 
induced by the tranverse field occur, the spontaneous magnetization M x is written in terms 
of the asymptotic behaviour of the determinant of a Toeplitz matrix which corresponds to 
the static two-spin correlation in the x direction, at large separation. The critical behaviour 
is determined, in high numerical accuracy, by adjusting the numerical results obtained for 
finite Toeplitz matrices by using a non-linear regression of the scaling relation and the 
e— algorithm. This is a rather remarkable result since it allows us to determine very precisely 
the critical behaviour of the system within a numerical approach. As expected, we have 
shown from these results that the inhomogeneous model belongs to the same universality 
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FIG. 13: (a) The real and (b) imaginary parts of the dynamic susceptibility Xq Z (^), at T = 0, for 
n = 4, m = /j,2 = Hz = /i4 = 1, Ji = 1, J2 = 1.5, J3 = 1.75, J4 = 2, h = 0.5 and 7 = 0.1, as 
function of frequency oj and for different values of q. 



class of the homogeneous one. 

For n greater than one, the system can present multiple phase transitions and it is shown 
that the divergence of the isothermal susceptibility x¥ > which is associated to inflection 
points in the induced magnetization, is also a signature of the quantum transitions. These 
critical points, as usual, correspond to the the points where the spontaneous magnetization, 
which is the order parameter characterizing the quantum transition, goes to zero. 

Explicit results are presented for n = 4 and n — 5, where we show that, differently from 
the isotropic model where we always have n transitions, the number of quantum transitions 
is dependent on the values of the anisotropy parameters 7. We have also concluded that, for 
n even, there are two critical 7's, whereas, for n odd, there just one critical 7 equal to zero, 
and this constitutes the main effect of the cell size on the critical quantum behaviour of the 
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anisotropic model. 

It has also been shown that the static two-spin correlation in the z direction, as in the 
homogeneous model, can present oscillatory or monotonic behaviour depending on the values 
of 7's and the field. The limiting surfaces, in the parameter space, which separates the two 
regimes, are the so-called disorder surfaces. For the special case when we have different 
J's and identical 7's, the disorder surfaces collapse into disorder lines and, in this case, 
they have been determine analytically. In particular, for n — 2, it has been proved exactly 
that these lines corresponds to regions where the model can be mapped onto an equivalent 
isotropic model. It should be noted that along the disorder lines the direct static correlation 
(rf (0)7^(0)) — (rf) 2 presents an exponential asymptotic behaviour, differently from the 
homogeneous model where it is equal to zero. 

We have also obtained the static and dynamic correlation on the field directions and, 
from these results, we have determined the dynamic wave-vector dependent susceptibility 
Xq Z {uj). Finally, we have shown that, as in the homogeneous model, the static suceptibility 
Xo 2 (0) is equal to the isothermal one \t ■ 
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IX. APPENDIX 

In this appendix we present analitycally the mapping of the anisotropic model onto the 
isotropic one, for a given set of parameters, for a chain with n — 2. This mapping will 



be obtained, as in the homogeneous model 25| . by imposing that the isotropic and the 



anisotropic models have the same spectrum. As it will be shown, this equivalence is only 
possible when 71 = 72. 

Therefore, let us consider the matrices A q and B 9 , given by eqs. ()14p and (j!5p . which in 
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this case are explicitly given by 



fiih f + fexp(-i2g) 



£ + f exp(i2g) 



fi 2 h 



(60) 



^-^exp(-i2g) 



■^ + ^exp(i2g) 







(61) 



From these results we obtain the matrix (A q + M q ) (A q — M q ) , whose eigenvalues, determined 
from eq. (|T7|) , which constitute the two branches of the spectrum of the chain which are given 
by 



1 1 

E q,k = l,2 = ^ l ai + " 2 + ( Cl + ° 2 ) COS 2 <?] ± 2 i ^ + " 2 + ( Cl + ° 2 ) COS 2< ?] Z 

- 4 [(ai + a cos 2g) (a 2 + c 2 cos 2g) - [3 2 - 5 2 - 2/35 cos 2q] } 1/2 



(62) 



where 



ol\ = [L-Ji + a x + b 2 ; a 2 = /i 2 /i + a 2 + &ij 
ai 



y(l+7i); a 2 = y(l + 7 2 ), 



6i = y (1 - 7i) 5 fc 2 = y (1 - 72) , 

d = M (1 + 7l ) (1 _ 72 ) ; C2 = M (1 _ (1 + 

(3 = + /i 2 ai)h; 5 = (/^a 2 + n 2 b 2 )h. 



(63) 



From the previous result we can obtain the spectrum of the isotropic chain, with parameters 
Hih , [i 2 h , J'i, J 2 , by making 71 = 7 2 = and from this we can write 



£g,jfc=l,2 



j'2 j'2 

(/if + fi 2 )h' 2 + -77- + -TT + j'xJ 2 cos2g 



± 



± l - +n' 2 )h' y/ (/4 " /4) 2 ^' 2 + + ^ 2 + 2 ^ J2 cos 2g. 



(64) 



By imposing the equivalence of the two spectra we conclude immediately that we should 
have ci = c 2 = c, which implies that 7l = 72 = 7. Introducing these conditions in eq. fl62|) 
we obtain 



^fe=i,2 = ^ ["l + "2 + 2c cos 2q] ± - ^ (e*i - a 2 ) 2 + 4 {(3 2 + 5 2 ) + 8/35 cos 2q, (65) 
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and by comparing it with eq. fl64|) we can write immediately the mapping relations 

(tf + rf)h' 2 + + f = (A*? + A^ 2 + (1 + 7 2 ) , (66) 

■AJ' 2 = JiJ 2 (1 - 7 2 ), (67) 

(a*? - // 2 2 ) 2 // 4 + // 2 (a4 + /4) 2 (j? + j; 2 ) = (a*? - ^)^ 4 + 

+7 2 (Jl - Jlf + (Jl + Jl) [(A*i + ^) 2 ^ - 7 2 (Mi - h 2 ] , (68) 
04 + M 2 ) 2 ^' 2 X4 = « [(Mi + /^ 2 ) 2 ^ - 7 2 (Mi " M 2 ) 2 h 2 } . (69) 

This set of equations is not independent, and this means that the parameters of the equiv- 
alent isotropic model are not uniquely determined. Moreover, we can show that eqs. (166]) 
to (1691) can only present a solution if the parameters of the anisotropic model satisfy the 
equation 

{Jl - Jlf (1 - 7 2 ) 2 - 8/iiM 2 (J? + Jl) (1 - 7 2 ) + 16A» 4 = 0. (70) 

This expression is obtained after some lengthy but straigthforward algebra, by using the 
result 

0*1 + a4) 2 h' 2 = ^ + ^?h 2 -i 2 ^-^?h^ 

1 — 7 2 

which is obtained from eqs. (l66l) and (|69l . and eqs. (!6"6l . (!6"7l) e ( JTB . The solution of eq.(|70l) 
gives the result 



\J x ±J 2 \Jl-j 



h = ± ! <*v = _^ = hc ^ 1 _ 2 (72 ) 
2 v //xi/x 2 

where /i c are the critical fields of the anisotropic model in the limit 7 — > 0. Although this 
result, for the disorder lines, is restricted to n = 2, we have verified numerically that, for 
identical 7's, it is still valid for arbitrary n. In particular, for n=4, the numerical verification 
of eq.(|72l) is shown in Fig.®. 
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7 


h c (exact) 


M 


h c (calculated) 


(3 






h c i 


= 0.130554559... 


0.479(8) 


0.13050(1) 


0.125(3) 


7 


= 0.1 




= 0.914163156... 


0.465(4) 


0.9142(1) 


0.125(2) 






h C 3 = 


= 1.310209705... 


0.436(8) 


1.31017(6) 


0.125(4) 






h C 4 


= 1.5922623078... 


0.395(2) 


1.59227(6) 


0.126(1) 






h c i 


= 0.9541874918... 


0.488(12) 


0.9542(3) 


0.125(6) 


7 


= 0.15 


h c2 


= 1.2687172332... 


0.465(8) 


1.26870(5) 


0.125(4) 






h C 3 


= 1.5905738964.. 


0.423(5) 


1.5906(1) 


0.125(3) 


7 


= 0.5 


h c -- 


= 1.5644253446... 


0.501(9) 


1.5644(2) 


0.125(4) 



TABLE I: Numerical results for the fitting, in the scaling region, of the spontaneous magnetization 
M x at T = 0, as a function of the uniform field for n = 4, [i\ = /j,2 = H3 = m = 1, 71 = 72 = 
73 = 74 = 7, J\ = 1, J2 = 1.5, J3 = 1.75, J4 = 2 and different 7. The parameters shown have been 
obtained from the non-linear regression of the equation M x (h) = Mq \h c — fo| 
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